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Abstract 

In the recent experiments designed to synthesize the element 115 in the 243 Am+ 48 Ca reaction at 

O 

Dubna in Russia, three similar decay chains consisting of five consecutive a-decays, and another 
different decay chain of four consecutive a-decays are detected, and the decay properties of these 
synthesized nuclei are claimed to be consistent with consecutive a-decays originating from the 
^vj ' parent isotopes of the new element 115, 288 115 and 287 115, respectively|l]. Here in the present work, 

the recently developed deformed RMF+BCS method with a density-independent delta-function 
interaction in the pairing channel is applied to the analysis of these newly synthesized superheavy 
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nuclei 288 115, 287 115, and their a-decay daughter nuclei. The calculated a-decay energies and 



m 

O 
O 

(N 



half-lives agree well with the experimental values and with those of the macroscopic-microscopic 
FRDM+FY and YPE+WS models. In the mean field Lagrangian, the TMA parameter set is used. 
Particular emphasis is paid on the influence to both the ground-state properties and energy surfaces 
introduced by different treatments of pairing. Two different effective interactions in the particle- 
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particle channel, i.e., the constant pairing and the density-independent delta- function interaction, 
together with the blocking effect are discussed in detail. 
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I. INTRODUCTION 

Since the prediction of the existence of superheavy islands in 1960s , the synthesis 
of superheavy elements has been a hot topic in nuclear physics. Following numerous ground 
breaking technical developments in the last three decades, the process of synthesizing 
superheavy elements has been sped up dramatically. From 1995 to 1996, Hofmann et al. 
DDDD at GSI in Germany successfully produced the elements Z=110, 111, and 112 by 
using low-energy heavy-ion collisions. In January 1999, the new element Z=114 was reported 
at Dubna in Russia 0,0]. Two years later, the element Z=116 was also reported at Dubna 
[loj | . In August 2003, in the reaction 243 Am+ 48 Ca held at Dubna Q], with a beam dose of 
4.3 x 10 18 248-MeV and 253-MeV 48 Ca projectiles, nine new odd-Z nuclei originating from 
the isotopes of the new element 115 , 288 115 and 287 115, were produced, respectively. So far, 
all elements with 110 < Z < 116 have been produced successfully in laboratory. All these 
exciting discoveries have greatly extended our knowledge about superheavy nuclei around 
the predicted superheavy islands and drawn more and more attention from the theoretical 
side. 

The experimental progress has led to a large-scale investigation of superheavy nuclei 
by both refined macroscopic-microscopic (MM) models such as the finite-range droplet 
model with folded- Yukawa single-particle potentials (FRDM+FY) ll( or the Yukawa-plus- 
exponential model with Woods-Saxon single-particle potentials (YPE+WS) . and micro- 
scopic models such as the Skvme- 
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tivistic mean field model[ 
element island", i.e, the next spherical doubly magic nucleus, is a dream for physicists for 
the past several decades. There are already many works in this respect (see Ref. 0, 14, ^ 
and references therein). Possible candidates predicted by different theories are 2 g 8 114, 2 72l20 
or even ^4126. However, due to the limit of proper projectiles, the superheavy elements syn- 
thesized are always neutron deficient and lie in the deformed region. The deformation effects 
are very hnportan, to understand the nuclear structures in superheavy nndei UUU- 
It is known experimentally that the heavy nuclei of the actinum series (Z=93-103) are well 
deformed and Bohr and Mottelson 0] also pointed out that the deformation can increase 
the stability of the heavy nuclei. The microscopic and self-consistent relativistic mean field 
model due to its natural description of spin-orbit interaction 0, 21,0], which is a purely 
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relativistic effect, has 
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3een proved to be a reliable method to describe exotic and superheavy 



In the present work, the recently developed deformed RMF+BCS method with a density- 
independent delta-function interaction in the pairing channel 23] is adopted to analyze 
properties of lately synthesized superheavy nuclei 288 115, 287 115 and their a-decay daughter 
nuclei. The density-independent (or dependent) delta- function interaction has been proved 
to be very successful to take into account the continuum effect both in relativistic and non- 
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28| . In the mean field part, 



is used, which has been proved to be very successful in describing 
18| | in the relativistic mean field model. Particular emphasis is 
put on the effective interactions used in the particle-particle channel and blocking effects. 
In what follows, we discuss in detail numerical details and results of our calculations. 



II. NUMERICAL DETAILS AND RESULTS 



The RMF calculations have been carried out using the model Lagrangian density with 



nonlinear terms both for the a and to mesons as described in detail in Ref . |23|, |29( , which 
is given by 
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where all symbols have their usual meanings. The corresponding Dirac equation for nucleons 
and Klein-Gordon equations for mesons obtained with the mean field approximation are 
solved by the expansion method on the widely-used axially-deformed Harmonic-Oscillator 



basis 23J, |30(. The number of basis used for expansion is chosen as Nf = Nf, = 20. More 
basis have been tested for convergence considerations. We use the parameter set, TMA {2^, 
for the RMF Lagrangian. 

Three kinds of approaches to take into account the pairing correlations have been adopted 
in the present work. The first is the usual RMF+BCS calculation with constant pairing 
interaction. The inputs of pairing gaps are A n = A p = 11.2/v^4 and the blocking effect 
is ignored. Such calculations have been performed extensively by Ren et al. [l8|. The 



second is the RMF+BCS calculation with a density-independent delta-function interaction, 
V = — Vo<5(r*i — F2) 23]. Here, the blocking effect is also ignored for comparison. The third 
is the same as the second one except that the blocking effect is taken into account by the 



usual blocking method 



23, 



31 



The pairing force strengths Vq are fixed by obtaining 
similar binding energy for 288 115 as the first approach, i.e., Vq = 280 MeV fm 3 in the second 
and Vq = 330 MeV fm 3 in the third calculations respectively. The same Vq has been used for 
both protons and neutrons. A slight change of the pairing strength, say ten percent, only 
changes the absolute binding energy less than 1.0 MeV and other results are hardly changed. 
Throughout the paper, the fist, second, and third kind of calculations will be abbreviated 
by "Const", "Deltal" and "Delta2". 



A. a-decay energies 

In Table. HJ a comparison for binding energies and a-decay energies between our 
three calculations, Const, Deltal, Delta2, the results from the macroscopic-microscopic 
FRDM+FY model Q and the experimental valnes for the -115 chain and the -'115 
chain is tabulated. More detailed properties obtained from the calculations Delta2 are 
shown in Table. |H] The theoretical half-lives T a are calculated with the well-known 
Viola-Seaborg formula [ll|,|33j]. The difference between the predicted Q a by Const, Deltal, 
Delta2, the FRDM+FY model Q, the YPE+WS model H and the experimental value, 
A a (theo.) = Q a (theo.) — Q a (expt.), is plotted in Fig. 1 and Fig. 2 for the 288 115 chain and 
the 287 115 chain , respectively. 

First, for the 288 115 chain, we notice that the quality of agreement between our calcula- 
tions ( particularly Delta2 ) and the experimental values is similar to those of MM models 
(FRDM+FY and YPE+WS). For the last two nuclei in the 288 115 chain, 272 1 07 and 276 1 09, 
results of MM models are closer to experimental values. For 280 111, our calculations are 
between the FRDM+FY model and the YPE+WS model. For 284 113, predicted a-decay 
energy by our calculations, similar to that of the YPE+WS model, is larger than the exper- 
imental value while the FRDM+FY model predicts a smaller value. The biggest difference, 
about 1.0 MeV, between theory and experiment is also found for this nucleus. For 288 115, 
both our calculations and the FRDM+FY model predict similar values that are smaller 
than the experimental value, while the result from the YPE+WS model is larger than the 
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experimental value. 

Second, for the 287 115 chain, similar things happen. For 271 107, Delta2 and the 
FRDM+FY model predict similar Q a , while the YPE+WS model predicts a larger value. 
Because no experimental value is observed for this nucleus, prediction of Delta2 is taken 
as the experimental value for comparison. For 275 109, predictions of all our three calcula- 
tions deviate from the experimental value more than those of the MM models. While for 
279 111 and 287 115, our calculations are closer to experimental values than the MM models. 
For 283 113, just like the case of 284 114, the difference between theory and experiment is 
relatively large. Our calculations and the YPE+WS model predict different trend for this 
nucleus from the FRDM+FY model also. 

Third, we note that all our three calculations predict similar a-decay energies for both 
the 288 115 chain and the 287 115 chain. The calculations Const and Deltal give essentially 
the same results for both decay chains. While the calculations Delta2 are generally better 
than the other two calculations. This is more obvious for the odd-odd 288 115 chain than 
for the odd-even 287 115 chain. Since the main difference between the second and the third 
calculations is the blocking effect, we could safely conclude that a proper blocking treatment 
can improve the calculated observables for odd-even or odd-odd nuclei. We also can say 
that the RMF+BCS calculations with constant pairing can describe even-even superheavy 
nuclei reasonably well. For odd-odd and odd-even nuclei, the inclusion of the blocking effect 
becomes necessary. We will see this point more clearly in the following discussions. 

Last, we would like to point out the difference between the MM models and the relativistic 
mean field theory used here. As we have seen in the above discussions (see also Table. |J 
and Fig. predictions of the MM models are closer to the experimental values for 

the elements 109 and 107, while for the other three elements, our calculations are better 
than those of the MM models. The reason could be that the MM models depend more on 
the knowledge of known nuclei. While the RMF model, due to its natural description of 
spin-orbit interaction, after including deformation, pairing interaction and blocking effect 
properly, could be more powerful in predicting the properties of unknown nuclei. 
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B. Energy surfaces and ground-state deformations 



Now let us discuss a bit more about the differences between our three different kinds 
of calculations. We have performed the constrained quadrupole calculations j^, 3^ for 
both the 288 115 chain and the 287 115 chain in all the three calculations. The corresponding 
energy curves are shown in Fig. 3 and Fig. 4. We should mention that such calculations are 
very time-consuming. First thing we see is that Const and Deltal give quite similar energy 
curves for both decay chains. In fact they also give almost the same results for all calculated 
quantities (see also Table. |U and Fig. 1-2), except for the a-decay energies where Deltal 
is better. Another noticeable difference is that the hight of the barrier between different 
minima can differ a little bit. In most cases, Deltal give shallower barriers than Const. 
Second we could see the difference between calculations with and without blocking, Delta2 
and Deltal, is relatively large, despite that the ground-state properties are quite similar. 
This shows that proper treatment of blocking effect is definitely necessary for odd-odd or 
odd-even nuclei (see also Fig. 0|2j). For the 287 115 chain, due to the way that we fixed 
the pairing strength Vo, the absolute binding energies from calculations with and without 
blocking differ around 1.0 MeV for some nuclei. 

Unlike medium or light nuclei where generally only two minima (one oblate minimum 
and one prolate minimum) or one spherical minimum are observed, the energy curves of 
superheavy nuclei are more complicated as we can see in Fig. 010 This is not surprising. 
As there are more levels in heavy nuclei, level crossing is more frequent to happen and 
lots of local minima may appear. Except for 288 115 and 284 113 in the 288 115 chain, 287 115 
and 283 113 in the 287 115 chain, the ground state of other nuclei can be determined without 
ambiguity, i.e., around (3 2 ~ 0.2. Similar results have been obtained by the macroscopic- 
microscopic YPE+WS model Q, more specifically, f3 2 =0.200, 0.211 and 0.224 for 280 111, 
276 1 09 and 272 107; (3 2 =0.202, 0.215 and 0.228 for 279 111, 275 1 09 and 271 107. The YPE+WS 
model predicts f3 2 = 0.138 and fl 2 = 0.149 for 284 113 and 283 113, which are also close to 
our calculations (3 2 = 0.17 and (3 2 = 0.18. While for 288 115 and 287 115, the YPE+WS 
model predicts [3 2 = 0.072 and (3 2 = 0.066 respectively, which are quite different from our 
calculations, (3 2 ~ 0.5 for both these two nuclei. This difference can be understood easily 
because these MM models predict 134 114 to be the next spherical doubly magic nucleus, while 
most self-consistent models shift this property to the more proton-rich side Q]. Further 
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experiments are needed to clarify these discrepancies between different theoretical models 
and different parameter sets in the same model. In our calculations, two other configurations 
(3 2 ~ —0.2 and /3 2 ~ 0.3 are also possible for 288 115 under certain conditions such as the case 
of a-decay. That is to say, decay from these two configurations to 283 113 are also possible 
based on the calculated a-decay energy. For 283 113, we find that the minima around /?2 ~ 0.2 
and @2 ~ 0.5 are close to each other. 

Since we see that isotopes of the element 115 are very deformed in our calculations, we 
would like to have a closer look at this element and the element 117, the mother element of 
the element 115 in the a-decay chain. The corresponding energy surfaces from all the three 
calculations are plotted in Fig. El for 292 117, 288 115, 291 117 and 287 115. It is clearly seen 
that the configuration around f3% ~ 0.5 is still stable against fission even for the element 
117, but the barrier is lowered greatly for the calculation Delta2 than the other two calcu- 
lations. Such an influence to the fission barrier introduced by the blocking effect has been 
demonstrated by Rutz et al. {3j| in the RMF model. Here, we notice that the adoption 
of the density-independent delta-function interaction instead of the constant pairing in the 
pairing channel further reduces the fission barrier. Further calculations by Delta2 show that 
a-decay energies of 292 117 and 291 117 are, respectively, 10.71 MeV (5=2076.62 MeV) and 
10.83 MeV (5=2053.36 MeV), with T a = 2.23 s and T a = 0.49 s. It would be very interest- 
ing to synthesize the nuclei 292 > 291 H7 and measure the a-decay chains, since our calculations 
predict that these nuclei would make a decays. 

To summarize, the constrained calculations show that the energy curves are relatively 
complicated and the predicted ground-state deformations by our calculations are close to 
those by the YPE+WS model. Further comparisons show that the fission barriers are quite 
different for our three calculations, especially for calculations with and without blocking. 
This suggests that to study superheavy nuclei more reliably one needs to use a more realistic 
effective interaction in the pairing channel other than the const pairing interaction and at 
the same time include the blocking effect properly. 

III. CONCLUSION 

With the recently developed deformed RMF+BCS method with a density-independent 
delta-function interaction in the pairing channel, properties of the lately synthesized su- 
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perheavy nuclei 288 115, 287 115 and their a-decay daughter nuclei in Dubna [l| have been 
studied. The TMA parameter set is used for the mean field Lagrangian. In the particle- 
particle channel, three different treatments are introduced, i.e., the usual constant pairing 
without blocking, the delta-function interaction with and without blocking. Constrained 
quadrupole calculations have been performed also for all these three calculations. Relatively 
complicated energy curves are observed for these superheavy nuclei, especially those nu- 
clei with proton number larger than 111. The calculated a-decay energies, Q a , are found to 
agree well with the experimentally observed values and also are close to those of macroscopic- 
microscopic FRDM+FY model and YPE+WS model. Predicted ground-state deformations 
agree well with those of macroscopic-microscopic YPE+WS model. For odd-odd and odd- 
even superheavy nuclei, which we have studied here, the inclusion of the blocking effect can 
improve the overall performance and is thought to be necessary based on our calculations. 
For purposes other than studying the ground-state properties, a more realistic interaction 
in the pairing channel, such as the density-independent delta-function interaction used in 
the present work, would be necessary. 
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TABLE I: The binding energies, B, and a-decay energies, Q a , of decay chains of 288 115 and 
287 115. Listed are the RMF+BCS calculations with constant pairing, Const, with the delta- 
function interaction without blocking, Deltal, and with Blocking, Delta2. FRDM+FY are results 
from the finite-range droplet model with folded Yukawa single-particle potentials nj. The last 
column is the experimental Q a from Dubna [J. All energies are in units of MeV. 



Nuclei 


Const 


Deltal 


Delta2 


FRDM+FY 


Exp er iment 




B 


Qa 


B 


Qa 


B 


Qa 


B 


Qa 


r) 

*°ia 


288 1 1 r 

llo 


2059.10 


9.78 


2058.80 


9.91 


2059.03 
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2059.12 


10.12 


10.61 ± U.U6 


284 1 i o 
llo 


2040.58 


11.22 


2040.41 


11.04 


2041.03 


10.74 


2040.95 


9.15 


10.15 ± 0.06 


280 1 i i 
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onoo crn 

2023.50 


in t n 
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-1 A /id 
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ono i o i 
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n qt _i_ n n^; 
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276 inn 

-L \J<J 


2005.70 


9.75 


2005.30 


9.73 


2005.66 


9.42 


2003.64 


9.93 


n ok _|_ n nc 


272 1 07 


1987.15 


8.16 


1986.73 


8.27 


1986.78 


8.60 


1985.27 


8.88 


9.15 ±0.06 


268 105 


1967.01 




1966.70 




1967.08 




1965.86 






287 115 


2051.88 


10.96 


2051.72 


10.82 


2053.36 


10.90 


2052.72 


10.25 


10.74 ±0.09 


283 113 


2034.54 


11.31 


2034.24 


11.19 


2035.96 


10.98 


2034.68 


9.35 


10.26 ±0.09 


279 in 


2017.55 


10.55 


2017.13 


10.55 


2018.64 


10.33 


2015.73 


10.92 


10.52 ±0.16 


275 1 09 


1999.80 


9.67 


1999.38 


9.67 


2000.67 


9.52 


1998.36 


10.06 


10.48 ± 0.09 


271 107 


1981.17 


8.18 


1980.75 


8.29 


1981.89 


8.65 


1980.13 


8.66 




267 105 


1961.05 




1960.74 




1962.24 




1960.49 







11 



TABLE II: The binding energies, B, neutron and proton quadrupole deformations, fo n and fo p , 
neutron and proton rms radii, R n and R p , the calculated alpha-decay energies and life-lives, Q a 
and T a , of superheavy nuclei on the alpha-decay chains of 288 115 and 287 115 from the calculations 
Delta2. The last two columns are experimental decay energies and liftimes. All energies are in 
units of MeV and all radii in units of Fermi. 



Nuclei 


B 


for, 


for, 


/ \ , , 


/ ) 

P 




-*- Oi 


Oi Jr v J 


TLfexDt) 

1 Oi V,^ Ir / 


288 115 


2059 03 


0.48 


0.50 


6.58 


6.41 


10.30 


6.86 s 


10 61 ± 06 


S7 + \f ms 

° — 30 


284 113 


2041 03 


0.17 


0.17 


6.37 


6.18 


10.74 


111 96 ms 


10 15 ± 06 


48 + £-? 8 , s 


280 m 


2023.47 


0.18 


0.19 


6.34 


6.15 


10.49 


118.97 ms 


9.87 ±0.06 


3.6+f| s 

— l.O 


276 109 


2005.66 


0.20 


0.20 


6.32 


6.12 


9.42 


25.08 s 


9.85 ± 0.06 


0.72«; 87 s 


272 107 


1986.78 


0.20 


0.21 


6.30 


6.09 


8.60 


1953.31 s 


9.15 ±0.06 




268 105 


1967.08 


0.21 


0.22 


6.28 


6.06 










287 115 


2053.36 


0.48 


0.50 


6.56 


6.41 


10.90 


80.57 ms 


10.74 ± 0.09 


ms 


283 113 


2035.96 


0.18 


0.18 


6.36 


6.18 


10.98 


12.76 ms 


10.26 ± 0.09 


1001^° ms 


279 m 


2018.64 


0.20 


0.20 


6.34 


6.15 


10.33 


142.07 ms 


10.52 ±0.16 


170^0° ms 


275 109 


2000.67 


0.21 


0.21 


6.32 


6.12 


9.52 


5.77 s 


10.48 ± 0.09 


9.7lf 4 s 


271 107 


1981.89 


0.21 


0.21 


6.29 


6.09 


8.65 


604.91 s 






267 105 


1962.24 


0.22 


0.22 


6.27 


6.06 
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FIG. 1: The difference between calculated Q Q (theo.) and experimental Q a (expt.), A a (theo.) = 
Q a (theo.) — Q a (expt.) } for the 288 115 chain is plotted against mass number A. 
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FIG. 2: The difference between calculated Q Q (theo.) and experimental Q Q (expt.), A a (theo.) = 
Q a (theo.) — Q a (expt.) } for the 287 115 chain is plotted against mass number A. 
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FIG. 3: The energy surfaces for 288 115, 284 113, 280 111, 276 109, 272 107 and 268 105 are plotted as 
a function of mass quadrupole deformation, for three calculations: Delta2 (solid line), Deltal 
(dashed line) and Const (dotted line). 
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FIG. 4: The energy surfaces for 287 115, 283 113, 279 111, 275 109, 271 107 and 267 105 are plotted as 
a function of mass quadrupole deformation, for three calculations: Delta2 (solid line), Deltal 
(dashed line) and Const (dotted line). 
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FIG. 5: The energy surfaces for 292 117, 288 115, 291 117, and 287 115 are plotted as a function of mass 
quadrupole deformation, for three calculations: Delta2 (solid line), Deltal (dashed line) and 
Const (dotted line). 
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